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Abstract In this paper a hybridized weak Galerkin (HWG) finite element method for solving the Stokes 
equations in the primary velocity-pressure formulation is introduced. The WG method uses weak functions 
and their weak derivatives which are defined as distributions. Weak functions and weak derivatives can be 
approximated by piecewise polynomials with various degrees. Different combination of polynomial spaces leads 
to different WG finite element methods, which makes WG methods highly flexible and efficient in practical 
computation. A Lagrange multiplier is introduced to provide a numerical approximation for certain derivatives 
of the exact solution. With this new feature, HWG method can be used to deal with jumps of the functions 
and their flux easily. Optimal order error estimate are established for the corresponding HWG finite element 
approximations for both primal variables and the Lagrange multiplier. A Schur complement formulation of the 
HWG method is derived for implementation purpose. The validity of the theoretical results is demonstrated in 
numerical tests. 
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1 Introduction 

Weak Galerkin (WG) refers to a general finite element technique for partial differential equations (PDEs) 
in which differential operators are approximated by their weak forms as distributions. Since their in¬ 
troduction, WG finite element methods have been applied successfully to the discretization of several 
classes of partial differential equations, e.g., second order elliptic equations the Bihar¬ 

monic equations [laiiaiiziisi], the Stokes equations [23], and the Brinkman equations [12]. WG method 
methods, by design, make use of discontinuous piecewise polynomials on finite element partitions with 
arbitrary shape of polygons and polyhedrons. Weak functions and weak derivatives can be approximated 
by piecewise polynomials with various degrees. The flexibility of WG method on these aspects of approx¬ 
imating polynomials makes it an excellent candidate for the numerical solution of incompressible flow 
problems. 

Hybridization of finite element methods is a technique where Lagrange multipliers are introduced to 
relax certain constrains such as some continuity requirements. The main feature of the HWG method is 
that their approximate solutions can be expressed in an element-by-element fashion. Hybridization [I] 
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can be employed to obtain an efficient implementation for solving PDEs. The generalization of this idea 
to mixed finite elements has been investigated in [2U51I20). The idea of hybridization was also used in 
discontinuous Galerkin methods [inilimS] to derive hybridizable discontinuous Galerkin (HDG) [7H5]. 

In this paper, the WG finite element formulation developed in [23] is hybridized to obtain our new 
hybridized weak Galerkin finite element method for solving Stokes equations. This HWG formulation 
can be modified easily to solve interface problems by adding two functionals arising from the jump 
condition to the right-hand side. A Schur complement formulation of the HWG method is derived for 
implementation purpose. By eliminating the interior unknowns and the Lagrange multipliers, the Schur 
complement formulation yields a system with much smaller size. We shall show that hybridization is a 
natural approach for the weak Galerkin finite element method of [53] . We shall also establish a theoretical 
foundation to address critical issues such as stability and convergence for the HWG finite element method. 

The paper is organized as follows. In Section 2, we briefly discuss the continuous Stokes problem 
and recall some basic results for later reference. After presenting some standard notations in Sobolev 
spaces in Section 3, we introduce two weakly-defined differential operators: weak gradient and weak 
divergence. The HWG finite element scheme for the Stokes problem is developed in Section 4. In Section 
5, we shall study the stability and solvability of the HWG scheme. In particular, the usual inf- sup 
condition is established for the HWG scheme. In Section 6, we shall derive an error equation for the 
HWG approximations. Optimal-order error estimates for the WG finite element approximations are also 
derived in this Section. The equivalence of HWG formulation and its Schur complement formulation is 
proved in Section 7. Finally in Section 8, numerical experiments are conducted. 


2 The model problem 

Let H C be a polygonal or polyhedral domain in for d = 2, 3 respectively. As a model for the 
flow of an incompressible viscous fluid confined in H, we consider the stationary Stokes problem with 


nonhomogeneous Dirichlet boudary conditions, given by 



— Au 4- Vp = 

f, 

in H, 

(2.1) 

V • u = 

0, 

in H, 

(2.2) 

u = 

g, 

on dfl. 

(2.3) 


Throughout the presentation, we assume that the unit external volumetric force acting on the fluid 

fe [l2(h)]‘^. 

The weak form in the primary velocity-pressure formulation for the Stokes problem (I2.1I) - (I2.3D seeks 
u € and p G Lo(^) satisfying u = g on dfl and 

(Vu,Vv)-(V-v,p) = (f,v), (2.4) 

(V-u,g) = 0, (2.5) 


for all V € [iLg(H)]‘^ and q G 

Recently a weak Galerkin finite element method has been developed for solving the Stokes equations 
in |23] . The main idea of weak Galerkin finite element methods is the introduction of weak functions 
and their corresponding weak derivatives in the algorithm design. With well defined weak functions and 
weak derivatives, a weak Galerkin finite element formulation for the Stokes equations is derived from 
the variational form of the PDF (I2.4I1 - (I2.5I1 by replacing regular derivatives with weak derivatives and 
possibly adding a parameter independent stabilizer: find Uh and ph from properly-defined finite element 
spaces satisfying 


(Vu,W, V^v) - (V™ • ^,Ph) + s{uh, v) 

(Vu, ■ Uh,q) 


(f>v), 

0 


( 2 . 6 ) 

(2.7) 
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for all test functions v and q in test spaces. In this paper, the WG finite element formulation developed 
in [23] is hybridized to obtain our new hybridized weak Galerkin finite element method for solving Stokes 
equation (I2.1I) - (I2.3I1 . 


3 Weak Differential Operators and Discrete Weak Gradient 


Let D be any open bounded domain with Lipschitz continuous boundary in = 2,3. We use the 

standard definition for the Sobolev space H^{D) and their associated inner products norms 

II • lls.D, and seminorms | • |s,d for any s ^ 0. For example, for any integer s ^ 0, the seminorm | • is 
given by 



with the usual notation 


d 

a = {ai,--- ,ad), \a\ = ai + -had, d°‘ = Y[dx-- 

i=i 

The Sobolev norm || • ||m,_D is given by 





The space H^{D) coincides with for which the norm and the inner product are denoted by 

II • ||d and (•, ■)d-, respectively. When D = fl, we shall drop the subscript D in the norm and in the inner 
product notation. 

The space H{div;D) is defined as the set of vector-valued functions on D which, together with their 
divergence, are square integrable; i.e.. 


iL(div; £») = {v : V e V • v e L'^{D)}. 


Let T be a polygonal or polyhedral domain with boundary dT. A weak vector-valued function on the 
region T refers to a vector-valued function v = {vo,Vb} such that Vq € [L^(T)]‘^ and Vf, € [Hi{dT)]‘^. 
Let 


V(T) = {v = {vo,v,} : Vo G [L^(r)]^v, G [H^dT)f} (3.1) 

Recall that, for any v G V(r), the weak gradient of v is defined as a linear functional Vu,v in the dual 
space of [H{div,T)]'^ whose action on each q G [H{div,T)]‘^ is given by 

(V^v, q)T = -(vo, V • g)T + (vfe, q ■ n)oT, (3.2) 

where n is the outer unit normal vector along dT, (•, ■)t stands for the L^-inner product in [L^(T)]^ and 
(•, ■)gT is the inner product in [Hi{dT)]‘^. 

A discrete version of the weak gradient operator is an approximation, denoted by in the 

space of polynomials of degree r such that 

(Vu,,r,TV,g)T =-(vo, V • q)t + (v6,g • n)aT, Vq € [PriT)f^‘^. (3.3) 

From the integration by parts, we have 


(vo, Vg)T = -(Vvo, q)T + (vo, q ■ n)oT 
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Substituting the above identity into (13.31) yields 

(Vu,,r,TV, q)T - (Vvo, q)T = (vb - Vo, g • (3.4) 

for all q € [Pr(T)Y^<^. 

To define a weak divergence, we require weak function v = {vq, v;,} to be such that vq € [L‘^{T)Y and 
Vb • n G H~^ {9T). Denote 

D(T) = {v = {vo, Vb} : Vo G [L\T)]f Vb • n G R-^dT)} (3.5) 

Recall that, for any v G V(r), the weak divergence of v is defined as a linear functional • v in the 
dual space of H^{T) whose action on each ip G H^{T) is given by 

(Vu, • V , p)t = -(vo, Vv?)t + (vb ■ n, p)dT, (3.6) 

where n is the outer unit normal vector along 9T, (•, ■)t stands for the L^-inner product in L'^(T) and 
{■,-)dT is the inner product in H^{dT). 

A discrete version of the weak divergence operator is an approximation, denoted by {yw,r,T‘) in 
the space of polynomials of degree r such that 

(Vu,,r,T ■ V, (/?)r =-(vo, V(/5 )t + (vb • n, (/7)aT, 'ip € PriT). (3.7) 

From the integration by parts, we have 

(Vo, V(y9)T = -(V • Vq, p)t + (Vq • n, p) dT ■ 

Substituting the above identity into (j3.3|) yields 

(V^,r,T • V , p)t - (V • \q,p)t = ((vb - Vo ) • n, p)qt (3.8) 

for all p G Pr{T). 


4 A Hybridized Weak Galerkin Formulation 

The goal of this section is to introduce a hybridized formulation for the weak Galerkin finite element 
algorithm that was first designed in [53]. 


4.1 Notations 


Let 7/t be a partition of the domain D into polygons in 2D or polyhedra in 3D. Assume that Th is shape 
regular in the sense as defined in [55]. Denote by £h the set of all edges or flat faces in Th, and let 
= £h\ d£l be the set of all interior edges or flat faces. Denote by hr the diameter oi T G Th and 
h = maxTeTh hr the meshsize for the partition Th- 

On each element T G Th, there are spaces of weak function V(r) and V{T) defined as in (13.11) and 
(13.51) . respectively. Denote by V and A the function space on Th and £h given respectively by 

A= n [i7^(5T)]". (4.1) 

TeTfe T&Th 

Note that the values of functions in the spaces V(ri) and V(T 2 ) are not related for any elements T and 
T 2 , even if T and T 2 share an interior edge or flat face e G £1°. The jump of v = {vq, Vbj on e is given by 


Me 


VblaTi - VblsTa, e G 
Vb, e G d£l. 


(4.2) 
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where VbldTi is the value of v on e as seen from the element Ti. The order of Ti and T 2 is non-essential 
as long as the difference is taken in a consistent way in all the formulas. Analogously, for any function 
A G A, we define its similarity on e G £h by 


((A))e 


A|ari+A|aT2, ^ G 

0, e G dfl. 


(4.3) 


Denote by ((A)) the similarity of X in £h- 

For any integer k ^ 1, denote by Wk(T) the discrete function space as follows: 

Wk(T) = {g:gG A^(D),g|T G Pk-i(T)}. 


Let 14 (T) denote the discrete weak function space as follows: 

Vk(T) = {v = {vo,v,} : {vo,Vb}|T G [Pk{T)f x [Pfc_i(e)]^e C dT}. 

Let Ak{dT) denote 

Afc(ar) = {A : A|e G [Pfe_i(e)]^e C OT}. 

By patching Wk{T), 14(T), and Ak{dT) over all the elements T G 7/t, we obtain three weak Galerkin 
finite element spaces Wh^ 14, and Ah given by 

Wh= X{ Wh{T), 14 = n = n ^k{dT). (4.4) 

TeTh TGTh TGTfc 

Denote by 14 ° the subspace of 14 consisting of discrete weak functions with vanishing boundary value 

14 ° = {v = {vo, Vb} G 14 : Vf, = 0 on 917}. 

Furthermore, let Vh be the subspace of 14 consisting of functions without jump on each interior edge or 
flat face 

Vb = {v G 14 : [v]e = 0, e G £1°}. 

Denote by Vj) a subspace of Vh consisting of functions with vanishing boundary values 

V° = {v G Vb : Vb|e = 0, e G 9D}. 

Let Sb be the subspace of Ah consisting of functions with similarity zero across each edge or flat face 


— {A G Ah : ((A))e — 0, e G £h}- 


The functions in the space S/i serve as Lagrange multipliers in hybridization methods. 

Denote by 'Vw,k-i and (Vu,_fc_i-) the discrete weak gradient and the discrete weak divergence on the 
finite element space 14. They can be computed by using (1531) and (13.71) on each element T, respectively. 

For each element T G Th, denote by Qq the projection operator from [L'^{T)]‘^ onto [Pfe(r)]'^. For 
each edge or face e G £h, denote by Qb the projection from [P^(e)]'^ onto [Pk-i{e)]‘^. We shall combine 
Qo with Qb by writing Qh = {Qo, Qb}- 


4.2 Algorithm 

On each element T G 74, we introduce four bilinear forms given below: 


st(v, w) 

= P4^(Q&vo - Vb,(3bVifo - Wb)aT, 

(4.5) 

aT(v, w) 

= (V^v, Vu,w)t + st(v, w). 

(4.6) 

&T(v,g) 

= (V^-v,g)T, 

(4.7) 

ct(v, A) 

= (vb,A)aT, 

(4.8) 
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for V = {vo, Vf,} G V'fc(r), w = {wo, Wf,} G VkiT), q G Wh{T) and A G Ak{dT). 
Their sums over all T G 7^ yield four globally-defined bilinear forms; 


s(v,w) 

= ^ st(v,w), v,wG14, 

T&Th 

(4.9) 

a(v, w) 

= ^ aT(v,w), v,wG14, 

TsTh 

(4.10) 

&(v,q) 

= bT{v,q), vGVh,qGWh, 

Tan 

(4.11) 

c(v. A) 

= ^ ct(v. A), vGVb,AGAb. 

Tan 

(4.12) 


The following weak Galerkin finite element scheme for the Stokes equation m was introduced and 
analyzed in [23] . 

Weak Galerkin Algorithm 1. A numerical approximation for (l2.1D - (l2.3p can be obtained by seeking 
u?i = {uq, Ub} G 14 and ph G Wh such that ut = Qfcg on dfl and 

a(ub,v) - 6 (v,pb) = (f,vo), (4.13) 

b{uh,q) = 0, (4.14) 


for all V = {vo, v;,} G and q G Wh- 

The weak Galerkin finite element algorithm 1 can be hybridized in the finite element space Vh by 
using a Lagrange multiplier that shall enforce the continuity of the functions in 14 on interior element 
boundaries. The corresponding formulation can be described as follows. 

Hybridized Weak Galerkin (HWG) Algorithm 1. A numerical approximation for (I2.1I) - (I2.3D can 
be obtained by seeking {uh]Ph', ^h) G 14 x Wh x Eh such that Uf, = Qbg on 511 and 

a{uh,v)-b{v,ph)-c{v,\h) = (f,vo), 

b{uh,q)+c{uh,tJ-) = 0 , 

for all V = {vo, V;,} G 14°, q G Wh and /r G Eh- 

Lemma 4.1. The WG finite element scheme has a unique solution. 

Proof. Let f = 0, we shall show that the solution of (I4.15|) - (I4.16I) is trivial. To this end, taking v = Uh, 
q = Ph^ and fi — \h and subtracting (14.161) from (j4.15l) we arrive at 


(4.15) 

(4.16) 


a{uh,Uh) = 0 . 

By the definition of a(-, •), we know Vu,Ub = 0 on each T G 7b, uq = u;, on each dT. 

By (13.81) and the fact that uj, = Uq on dT we have, for any r G [Pk-i{T)]^^^, 

0 = (Vu,Ub,T)T 

= (Vuo,t)t - (uo - ut„r • n)aT 
= (Vuo,t)t, 

which implies Vuq = 0 on each T G Th and thus Uq is a constant. Since Uq = Uf, on each dT, we have 
b{uh, q) = -(uo, Vq) -b (uq • n, q) = (V • Uh, q) = 0. 


From (14.161) . we obtain 


c(ub,At) = 0. 
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Let jj, = |ufe] in above equation, then it follows that 

E /^IIMIIe=0. 

Ten 

Thus |uf,] = 0, which implies that Uq is continuous and we arrive at = {0,0} in il. 

Let Vf, = 0 in (I4.15L then it follows from ujj = (Q, Qj and f = Q that 

b{v,Ph) = (Vu, • v,Ph) = -(vo,Vph) = 0. 

Hence we have Vph = 0 on each T € Th- Thus ph is a constant in T. Let Vf,|e = |p?t]e)Vo = 0. Then 

c(v,Ah) = 0. 

Thus 

0 = b(v,Ph) = E WlPhlWl- 

ee£h 

Hence ph is continuous. From ph G Lg(H), we would obtain = 0 in H. 

Finally, Let = Xh, from u/j = {Q, Q} and = 0 in H, we obtain 

c{v,Xh) = 0, 

which means A/i = 0. 

This completes the proof of the lemma. □ 

4.3 The Relation between WG and HWG 


The rest of this section will show that the above two schemes are equivalent in that the solutions VL^^ph 
from (I4.13|) - (l4.14p and Ufi,Ph from (I4.15I) - (I4.16I) are identical, respectively. 

For any v € V°, let 

|||v|||^ = a(v,v) = E l|Vu,v|||+ E - Vfelllj.. (4.17) 

Ten Ten 

It has been verified in that (14.171) defines a norm in the vector space V°. 

Theorem 4.2. Let u/j £ Vh, pn G Wh be the first two components of the solution of the hybridized WG 
algorithm Then, we have |uje = 0 for all e £ i.e., u £ V/j and uj, = Qtg on dfl. 

Furthermore, we have that and ph satisfy the equation that is, u/j = u;^ and ph = Ph- 

Proof. Let e £ fj) be an interior edge shared by Ti and T 2 . By letting q = 0, fi = |u/i]e on e £ dTi, 
pb = — |u/j]e for e £ dT 2 , and = 0 otherwise in (I4.16L we have from (14.121) that 

0 = c{Uh,fJ,) = E {^h,fJ-)aT = / luhfeds, 

Ten 

which implies [u/^Je = 0 for any e £ 

Next by letting pt = 0, we obtain from (14.121) that 

b{uh,q) = 0 


for all q £ Wh- 

For any v £ V){, it follows from |v]e = 0 on any e £ fj) and v = 0 on dft that 


Thus, we arrive at 


C(V, Xh) = E i^b, K)dT = E (M’ = 0- 

Ten Ten 

a(u;^,v) - b{-v,ph) = (f,vo). 


which is the same as (I4.13p . It implies that {uh;ph) is a solution of the WG scheme (I4.13I) - (I4.14I) . It 
follows from the uniqueness of solution of (I4.13p - ()4.14p that u/i = and ph = Ph, which completes the 
proof. □ 
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5 Stability Conditions for HWG 

It is easy to see that the following defines a norm in the finite element space 


|A||h. = ^ he\\X\\ 


(5.1) 


As to V°, for any v let 


|v|lyo = 


E 


(5.2) 




We claim that || • \\yo defines a norm in V^. In fact, if ||v||yo =0, then |v]e = 0 on each interior edge or 
flat face e € £1°, and hence v £ Vj). Since ||| • ||| defines a norm in the vector space Vj), then v = 0. This 
verifies the positivity property of || • ||yo. The other properties for a norm can be checked trivially. 

Remark 5.1. Similarly, for any (p = (p; A) £ Mh, we can define 


UIIm, = Ibll 


lAlh,. 


(5.3) 


Lemma 5.2. ( [22) 1 (Trace Inequality) Let Th he a partition of the domain SI into polygons in 2D or 
polyhedra in 3D. Assume that the partition 7h satisfies the assumptions Al, A2, and AS as specified 
in J22^ . Let p > 1 be any real number. Then, there exists a eonstant C such that for any T £ Th and 
edge/face e £ dT, we have 

lblli.(e) ^ + /.^||V0|li.(^)), (5.4) 

where 9 £ IT^’^(T) is any function. 

Lemma 5.3. ( [22) 1 (Inverse Inequality) LetTh be a partition of the domain Ll into polygons or polyhedra. 
Assume that Th satisfies all the assumptions A1-A4 and p ^ 1 be any real number. Then, there exists a 
constant C(n) such that 


IIVpIIt.p ^ C(n)hf^y\\T,p, VT £ Th 

for any piecewise polynomial ip of degree n on Th- 

Lemma 5.4. (Boundedness) There exists a constant C > 0 such that 


(5.5) 


a(w,v)| 


C'||'w|bo||v|bo, 

V’w,v £ 14 °, 

(5.6) 

l&(v,9)l 

< 

C||v|bo|b||, 

yvevlqeWh, 

(5.7) 

|c(v,A)| 

< 

C||v|bo||A||H„ 

Vv £ , X € ^h- 

(5.8) 


Proof. To derive (j5.6|) . we use the Cauchy-Schwarz inequality to obtain 

|a(w,v)| = ^ (Vu,w, V^v)t + hf~^{Qtwo - '^b,Qb^o - ^b)dT 
TdTh 


XTeTfe / 


KTdTh 1 


+ E - Wfcl 

\T€Th 

^C'||w|lyo||v||yO. 


dT 


^ ^ hrp llQfcVO 
yTeTh 


'VbWdT 
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As to (Ell) , we use (|3.6L trace inequality (15.41) . and inverse inequality (15.51) to obtain 




^ (V^ • v,q)q 

TdTh 

- (’^o,Vp)y+ ^ {vt,pn)9T 

Ten Ten 

X] (^■vo,4')t- {^0-^b,pn)oT 


Ten 


Ten 


< (Eii^-^oIItJ (Xmi 

XTen / yTen 


+ ( X l|vo-Vbllarj [ X Ill’ll 

XTen ) yren 


9T 




C X Ml 

xren / 


nCh 2 I ^ llvo-VbllaTl (Ibll+/i||Vp||) 
\Ten / 

^ M|v|lyo|H|. 


As to (15.81) . it follows from the Cauchy-Schwarz inequality that 


|c(v,A)| = 


X dT 

Ten 


< 


X 


'IIM ,"2 




X 


^M|v|bo||A||H„ 


which completes the proof. 

Lemma 5.5. (Coercivity) For any v G V°, we have 

|a(v,v)| > C||v||^o. 


□ 


(5.9) 


Proof. For any v G V°, we have ||v||^o = |||v|||, which means the estimate (15.9p holds true. This completes 
the proof. □ 


Lemma 5.6. ( [23 ] ) There exists a positive constant j3 independent of h such that 

&(v,p) 


sup 
vev“ lll''^ll 




(5.10) 


for all p G Wh- 

Lemma 5.7. For any given p G Wh, there exist a positive constant fd independent of h and a v G 
such that 


Mllv^o 




(5.11) 
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Proof. From C and Lemma 15.61 we have, for any p G Wh , there exists a v S 


H^,P) ^ b{^, P) 
l|v||vO ll|v||| 




(5.12) 


which completes the proof of the lemma. □ 

Lemma 5.8. For any given t G there exist a v G with Vq = d and a constant C > 0 such that 

c(v,t) 


|v||vo 




(5.13) 


Proof. For any r G Eh, we have (('r))e = 0 or equivalently = 0 on each interior edge e £ £1° and 

X = 0 on all boundary edges. By letting v = {0, hf-r} G Vh in c(v, r) and s(v, v), we obtain 


and 


H + (v2,t 2), = 2 ^ hellxll^, 


s(v,v) =2 hj}hl{\\T'^\\l + ||t2||2) < 2 ^ he\\T\\l 


(5.14) 


(5.15) 


It follows from (13.31) . Cauchy-Schwarz inequality, the trace inequality (15.41) . and the inverse inequality 


that 


(V„v,V^v)t = ^ (v:,V^v)e s: X! ^e||T*|U||V^v||, 

e^dT eedT 

Y /ll||r*||e||V^v||T, 

egOT 


(5.16) 


where is chosen to be or according to the relative position of and e, which implies that 

lie. (5.17) 


IV^vIIt GiCY^i 

egar 


By summing over all elements, we obtain 


(V^V,V^V);, /le||T*||i 

It follows from (I5.15P and (15.181) that 

lllvf he||T*||2 = C||r|||^. 

By combining (15.141) and (15.191) . we obtain that there exists a constant C > 0 such that 

c(v,x) 


|v|lyo 




which completes the proof. 


(5.18) 

(5.19) 

(5.20) 

□ 


6 Error Estimates 

The goal of this section is to derive an error equation for the HWG finite element solution obtained from 
(I4.15p - (I4.16|) . This error equation is critical in convergence analysis. 

In addition to the projection Qh = {Qo, Qb} defined in the previous section, let Qh and Qh be two 
local projections onto Pk-i{T) and [Pk-i{T)]^^'^, respectively. 
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Lemma 6.1. ( [23] ^ The projection operators Qh, Q/i, and Qh satisfy the following eommutative prop¬ 
erties 


^n,{Qh^) = Q/.(Vv), Vve (6.1) 

• (Q/^v) = Q/i(V-v), Vv e iJ(div, 11). (6.2) 

Denote by (u;p) the exact solution of (I2.1I1 - (I2.3I1 . Let {\ih]Ph \ ^h) £ 14 x Wh x be the solutions of 
(I4.15|) - (l4.16p . Let A = Vu • n — pn. Define error functions as follows 

e/t = {Qou-uo,(5bU-Ufc}, eh = Qhp-ph, 8h = Qb\-\h- (6.3) 

Lemma 6.2. Let (u;p) he the exact solution of 112. and {uh',Ph', ^h) £ 14 x Wh x Eh be the 
solutions of Then, the error functions eh, Eh, and 6h satisfy the following equations 

a{eh,v) - b{-v,eh) - c{-v,dh) = ^u,p(v), V\e€Vh, (6.4) 

b{eh,q) -\-c{eh,fJ.) = 0, \/q € Wh, fx € Eh, (6.5) 


where 

4,p(v) = X] “ Qb(Vu)) • n)aT 

TdTh 

- X] ^'"0 - Vfe, {p - Q?iP)n)aT + s(QhU, v). 
TdTh 

Proof. First, applying (13.2L Lemma l6.11 and the integration by parts, we have 


w (^h^) , ^w^)t 


Summing over all T £ 7b reaches 


(Qb(Vu), Vu,v)t 

-(vo, V • Qh(Vu))T + (vb, Qh(Vu) ■ n)gT 
(Vvo,Qh(Vu))T - (vo - Vb, Qb(Vu) • njoT 
(Vvo, Vu)t - (vo - Vb, Qb(Vu) • njoT 
-(Au, Vo)t + (vo - Vb, (Vu - Qb(Vu)) • n)gT 
-l-(vb,Vu ■ n)gT- 


-(Au,vo) = (V^(Qhu),Viav) - (vo-Vb,(Vu-Qh(Vu))-njoT 

Ten 

~ X ('^b,Vu ■ n)oT- 
Ten 


Similarly, by using (13.61) and the integration by parts, we have 


(V™ • v,Qbp)T 


-(vo, V(Qbp))T + {vb, (Qbp)n)aT 

(V • Vo, QbP)T - (vo - Vb, (QhP)n)gT 

(V ■ vo,p)t - (vq - Vb, (QhP)n)gT 

-(vo, Vp)t + (vo,pn}gT - (vo - Vb, (QhP)n)gT 

-(vo, Vp)t + (vo - Vb, {p - QhP)n)gT + {vb, Vu • n)gT, 


( 6 . 6 ) 


(6.7) 


= {^w ■ V, QhP) + X ^'"0 - Vh, {p - QhP)n)gT 

Ten 

+ X ^'^b,prf}gT- 
Ten 


Summing over all T leads to 
(Vp,vo) 


( 6 . 8 ) 
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By using the identity —(Au, vq) + (Vp, vq) = (f,vo) and noticing that 

^ (vf,, Vu • n - pn)QT = c(v, A), 

TdTh 

we obtain 

a(v, Q;,u) - 6(v, Q_hP) - c(v, A) = (f, vq) + ^u,p(v). (6.9) 

Combining with the scheme (14.151) as follows 

a(v, u?,) - b{x,ph) - c(v, \h) = (f, Vo), 

we obtain 


a{eh,'v) - b{-v,eh) - c{\r,Sh) = ^u.p(v). 
As to (16.51) . from Theorem 14.21 we know that |e/j] = 0, which leads to 

c(e/,,/x) = 0, V/xGS^. 


Moreover, for any q G Wh, we have 

b{<^h,q) 


b{QhU, q) 

(V„ • (Q;^u),g)T 

T&Th 

Y {Qhi^ ■ u),q)T 

Ten 

(V-u,g) =0. 


This completes the proof. □ 

Next we shall establish some error estimates for the hybridized WG finite element solution {uh;ph', Ah) 
arising from (14.151) - (14.161) . The error equations ()6.4|) - (j6.5|) imply 

a{QhU-Uh,v) - b{v,QhP-Ph) - c{v,QbX-\h) = ■^u.p(v), Vv G V”, 

b{QhU-Uh,q) + c{QhU-Uh,r) = 0 , Wq £ Wh, pt &‘^h, 

where £u,p(v) is given by (1^^ . The above is a saddle point problem for which the Brezzi’s theorem [3] 
can be applied for an analysis on its stability and solvability. Note that all the conditions of Brezzis 
theorem have been verified in Section [S] (see Lemmas 15.4115.51 and 15.81) . The following error estimate can 
be proved similarly with Theorem 7.1 of |23j . 

Theorem 6.3. Let (u;p) be the exact solution of Ii2.1\) - li2.3\) and {uh',Ph',Xh) G 14 x Wh x Eh be the 
solutions of ^.15\ )- fi4.16] l. Then, there exists a constant C such that 

IIQhU- Uhliyo + \\QhP-Ph\\ + lIQbA- Ahlls^ 4 C'h'=(||u||fe+i + ||p||fc). (6.10) 

Theorem 6.4. Let (u;p) be the exact solution of S2. and Ah G Eh be the last component of the 
solution of On the set of interior edges £^, let A = Vu • n — pn. Then, there exists a 

constant C such that 

\\\-\hh,^Ch\\\uU+i + \\ph). 

Proof. From the triangle inequality and Theorem 16.31 we have 

IIA-Ahlb. 4 IIA-QhAlls. + IIQhA-Ahb., 

IIQhA-AhllH, 4 C'h'=(||u||h+i + |blU). 


( 6 . 11 ) 
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Thus we just need to concentrate on ||A — QbA||H;.- 

Applying Definition 15.11 trace inequality, and the property of LF' projection, yields 


IIA-QbAIII, = 




K\\\-Q,X\\l 

^ /ill Vu - Q,Vu||2 + ^ h\\p - <^up\\l 

ee£» ee£« 

C Y, (IIVu - Q,,Vu||2, + /i2||Vu - Qi,Vu||?,r) 
Ten 


+C' 'Y^ (lb ~ QhpWr + ^^Ib ~ QiipIIi.t) 
Ten 

< Ch^\\\uU+i + \\p\\k)\ 


which completes the proof. □ 

The following L^-error estimate for Qqu — Uq follows from Theorem 14.21 and Theorem 7.2 of p5] . 

Theorem 6.5. ( [53]) Let (u;p) with fc ^ 1 and {uh;ph',Xh) G 14 x 114 x S/j be the exact solution of 
M-ifO) and be the solutions of |^.75[ )- |4/d[ ), respectively. Then, the following optimal order error 
estimate holds true 

IIQou - uqII 4 C/i'=+i(||u||fc+i + llpllfe). (6.12) 


7 Efficient Implementation via Variable Reduction 

The degrees of freedom in the WG algorithm (I4.13|) - (I4.14|) can be divided into two classes: (1) the interior 
variables representing Uq, and (2) the interface variables for Uf,. For the HWG algorithm (14.151) - (14.161) . 
more unknowns must be added to the picture from the Lagrange multiplier A/j. Thus, the size of the 
discrete system arising from either (|4.13l) - (|4.14l) or ()4.15p - ()4.16|l is enormously large. 

The goal of this section is to present a Schur complement formulation for the WG algorithm (14.131) - 
(I4.14[) based on the hybridized formulation (|4.15I) - (I4.16I) . The method shall eliminate all the interior 
unknowns associated with Uq and the interface unknow \h, and produce a much reduced system of linear 
equations involving only the unknowns representing the interface variables Uf,. 

7.1 Theory of variable reduction 

Denote by Bh the interface finite element space defined as the restriction of Vh on the set of edges £h', 
i.e., 

= {v = {pL]p} : H G [Pfc_i(e)]'^,p|e G Pk-i{e),e G £h}. 

Bh is a Hilbert space with the following inner product 

{'^b,q.b)£h = X! Vw{,,qb G Bh- 

ee£h 

Let 5° be the subspace of Bh consisting of functions with vanishing boundary value on d£t. It is not hard 
to see that the interface finite element space Bh is isomorphic to the space of Lagrange multiplier 
The Schur complement through an elimination of the Lagrange multiplier Xh and the interior unknown 
Uq can be implemented through a map, denoted by 5'f. 

We define the map Sf : Bh ^ B^ as follows: for a fixed ph and any given function Wf, G Bh, the image 
£fi'^b',Ph) can be obtained by 

Step 1. On each element T G Tk, solve for wq in term of and ph from the following local problem: 

aT(w,i, v) - 5 t(v,p,,) = (f, vo)t, Vv = {vq, 0} G yfc(r), (7.1) 

where Wh = {wo,W(,} G Vk{T),ph G Wk{T). Denote -wq = Df{'Wb;Ph)- 
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step 2. On each element T € Th, solve for Ch,T & Ak{dT) in term of = {wojWb} and ph from the 
following local problem: 

ct(v, C/t,T) = aT(w/,,v) - Vv = {0, v4 £ 14(T), (7.2) 

Thus we obtain a function Ch,T S A/j. Denote Ch,T = 

Step 3. Define >S'f(w{,;p^) by the similarity of Ch on interior edges and zero on boundary edges; i.e., 

S{{wt]Ph) = iiCh))- (7.3) 

It follows from (17.31) that Sf{wi,]Ph) & The following are two properties regarding the operator Sf 
and the related terms: 

(1) The sum of (17.ip and (17.21) yields 

ct(v, C/i.t) = aT(w/,,v) - 6t(v,_p/i) - (f,Vo)T, Vv = {vq, v;,} G Vfc(r). (7.4) 


(2) It follows from the superposition principle we have the following identify 


-S'f(wb;p?i) = 5'Q(wb;p?,) + S'f(Q;0), \/wb G Bh,Ph & Wh, (7.5) 


where Sq is the operator with respect to f = 0. 

Lemma 7.1. The following identity holds true for the operator Sq: 


{SQ{wi,;ph),(li,)s^ = a{wh,qh) - b{(ih,ph), Vwf„qf,GB°, 


(7.6) 


where Wh = {Dg{wb-,ph),Wb} and q/^ = {DPh), <lb}■ 

Proof. For any wj,, q?, G B^, from the definition of the operator Sf we obtain 

w/t = {DQ{wb;ph),'v^b}, Ch =LQ{wb;Ph), 

q/i = {Do(q&;p/i),qb}. 

Letting f = 0 in (17.41) yields 


{SQ{wb;ph),<ib)eH 


E iCh,T,(ib) dT 
Ten 

E CT(q/i, C?i,t) 

Ten 

E -bT{qh,Ph)- 

Ten 


This completes the identity (17.61) . □ 

Lemma 7.2. LetUh = {uo,Ub} G Vh, Ph S Wh, and \h G 'E.h be the unique solution of the hybridized WG 
algorithm Then, we have G Vh and u;, G B^ is a well defined function. Furthermore, it 

satisfies the following equation 

Sf{ub-,Ph) = iiCh)) = Q. (7.7) 

Proof. Since {uh;ph',A.h) is the unique solution of the HWG scheme (I4.15I) - (I4.16L then from Theorem 
14.21 we have |u^]e = 0 on each interior edge or flat face e G f° and Uf, = Qbg on dfl. Thus, u/j G Vh and 
its restriction on £h is a well-defined function in Bh. 
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In order to verify (ITTI) . we take v = {vp, 0} £ Vk(T) on T and zero elsewhere in (14.1511 . it follows that 
u/j satisfies the local equation 


aT{uh,v) - bT{v,ph) = (f,vo)T, Vv = {vo,Q} £ Vk{T). 

Next, taking v = {Q,Vf,} £ 14(T) on T and zero elsewhere in (I4.15L yields that Xh satisfies the local 
equation 

CTiv,Xh,T) = aTi'^h,^) - bT{v,ph), Vv = {Q,vt,} £ Vfc(r), 
where Xh^T is the restriction of X^ on dT. Thus, from the definition of the operator Sf, we obtain 


Sf{ut-,ph) = {{Xh)). 


Combining with the fact that Xh £ S/j, we have {{Ck)) = Q, which completes the proof of the lemma. □ 

Lemma 7.3. Let £ Bh he a function satisfying Ub = QbO on dLl, u;, and ph satisfy the following 
operator equation 

Sf{xih-,Ph) = Q. (7.8) 

Then, {uh]Ph) £ 14 x Wh is the solution of the WG finite element solution arising from 
Here Up is the solution of the following local problems on each element T G 7h, 

aT{uh,v) - bT{v,ph) = if,vo)T, Vv = {vp, 0} £ 14(T), (7.9) 

with u/i = {up, Ub}. 

Proof. For each T GTh, solve for Xh^r £ ^k{dT) from the local problem 

CT{v,Xh,T) = aT{uh,v) - bT{v,ph), Vv = {0, Vb} £ Vfc(r). (7.10) 


Define a function Xh £ Xh by Xh\dT = Xh^r. Since (ub;p/t) £ Bh x Wh satishes the operator equation 
(I7.8[l . Ub satisfies the boundary condition, and Up is given by (I7.9|l . it follows from the definition of the 
operator Sf that 

{{Xh)) = Sf{uh;ph)=0, (7.11) 

which implies Xh £ 'E.h. 

By subtracting (17.101) from (17.9L we have 


aT(Ub,v) - bT{v,Ph) - CT{v,Xh,T) = (f,Vp)T, Vv = {vo,Vb} £ 14(T). 
By summing up the above equations over all T £ 1 ) 1 , we obtain 

a(ub,v) - b{v,ph) - c(v, Ab) = (f,vo), Vv = {vp,Vb} £ 44- 
By restricting v to the weak finite element space and using (17.1111 we arrive at 

C ( V , Xh) = ^ {{{Xh))e, Vb)e = 0. 


(7.12) 


(7.13) 


Thus we obtain 

a(ub,v) - b{v,ph) = (f,vp), Vv = {vo,Vb} £ 14°- 

Recall that the assumption U 5 = Qbg on dfl and Theorem 14.21 we have Ub is the WG finite element 
solution of (j4.13ll - (j4.14ll . This completes the proof of the lemma. □ 

Combining the above two lemmas yields the following result. 

Theorem 7.4. Let Ub £ Bh be any function such that Ub = QbS on dQ. Let Up be the solution of ([7.ff[ ). 
Then {uh',Ph) 'is the solution of if and only if Ub satisfies the following operator equation 


5'f(ub;pb) = Q. 


(7.14) 
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7.2 Computational algorithm with reduced variables 

Together with the equation (17.51) . (17.141) gives rise to 

S'g(ufc;p^) = -^f(Q;0). (7.15) 

Let Gf, G Bh be a finite element function such that Gf, = Qbg on the boundary of VL and zero elsewhere. 
From the linearity of operator S'q, we have 

Ph) ^biPh) 4“ Ph) ■ 

Substituting this equation into (17.151) . one obtains 

<S'Q(uft — Gb'yPh) = —-Sf (0; 0) — S'g(Gh; ph). 

Note that the function H;, = ut — Gt has vanishing boundary value. Letting rj, = — S'f (0; 0) — SQ{Gb]Ph), 
we have 

SQillb-,ph) = rb. (7.16) 

The reduced system of linear equations (17.161) is actually a Schur complement formulation for the WG 
finite element scheme (I4.13|) - (l4.14p . Note that (|7.16l) involves only the variables representing the value of 
the function on fj). This is clearly a significant reduction on the size of the linear system that has to be 
solved in the WG finite element method. 

Variable Reduction Algorithm 1. The solution {ub',Ph) to the WG algorithm (14.131) - (14.141) can be 
obtained in the following steps: 

Step 1. On each element T G Th, solve for r^ from the following equation: 

ff) = —-S'f(Q;0) — SQ{Gb',Ph) 

This step requires the inversion of local stiffness matrices and can be accomplished in parallel. The 
computational complexity is linear with respect to the number of unknowns. 

Step 2. Solve for by means of the operator equation (17.161) . 

Step 3. Gompute Uf, = G;, + H;, to get the solution on element boundaries. Then on each element T, 
compute Uq = Df{ub;ph) by solving the local problem (17.11) . This task can also be implemented in 
parallel, and the computational complexity is proportional to the number of unknowns. 

Note that, Step (2) in the Variable Reduction Algorithm 1 is the only computation extensive part of 
the implementation. 

8 Numerical Experiments 

The goal of this section is to report some numerical results for the hybridization weak Galerkin finite 
element method proposed and analyzed in previous sections. 

A Schur complement technique of the HWG method is utilized to decrease the degree of freedom. For 
example, if 11 = [0,1] x [0,1] and the uniform triangulation is used with mesh size V^/n, the number 
of elements is Nt = 2n^ and the number of edges is Ne = 3n^ + 2n. If fc = 1, then the degree of 
freedom for usual weak Galerkin method is 7Nt + 2Ne = 20n^ -b 4n, the degree of freedom for hybridized 
weak Galerkin method is 13Nt + 2Ne = 32n^ + 4n, while the degree of freedom can be reduced to 
2Ne + Nt = 8n^ + 4n by using the Schur complement. 

Let (u;p) be the exact solution of (I2.1D - (I2.3I) and {uh\Ph) be the numerical solution of (14.131) - (14.141) . 
Denote eh = and Eh = QhP — Ph- The error for the weak Galerkin solution is measured in four 

norms defined as follows: 

ll|e?i||P = X] ( / |Vu,e/ipdT-b/ (eo - eh)^dsV 

X" ^dT / 
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Example 7.1 Consider the problem (j2.ip - (l2.3|) in the square domain 17 = (0,1)^. The HWG finite 
element space fc = 1 is employed in the numerical discretization. It has the analytic solution 

C sin(27rx) cos(27r?/) \ 2 2 1 

and p = X y — -. 

— cos(27rx) sin(27r?/) J 9 

The right hand side function f in (12.11) is computed to match the exact solution. The mesh size is denoted 
by h. 

Table 7.1 shows that the errors and convergence rates of Example 7.1 in ||| • |||— norm and norm for 
the HWG-FEM solution u are of order 0{h) and 0{h?) when fc = 1, respectively. 

Table 7.2 shows that the errors and orders of Example 7.1 in norm for pressure and A . The 
numerical results are also consistent with theory for these two cases. 


Table 7.1. Numerical errors and orders for u of Example 7.1. 


h 

llle.i 

order 

lle.ll 

order 

1/4 

5.8950e+00 


1.3555e+00 


1/8 

2.9253e+00 

1.0109 

2.3750e-01 

2.5128 

1/16 

1.4552e+00 

1.0074 

4.9049e-02 

2.2756 

1/32 

7.2651e-01 

1.0022 

1.1500e-02 

2.0926 

1/64 

3.6312e-01 

1.0006 

2.8254e-03 

2.0251 

1/128 

1.8154e-01 

1.0001 

7.0325e-04 

2.0063 


Table 7.2. Numerical errors and orders for p and A of Example 7.1. 


h 

lk/.ll 

order 

IIQfeA-A^II 

order 

1/4 

5.1609e-01 


8.6908e-01 


1/8 

2.9426e-01 

0.8105 

3.2951e-01 

1.3992 

1/16 

1.4706e-01 

1.0007 

9.8958e-02 

1.7354 

1/32 

7.2990e-02 

1.0107 

2.6698e-02 

1.8901 

1/64 

3.6391e-02 

1.0041 

6.9159e-03 

1.9487 

1/128 

1.8180e-02 

1.0012 

1.7681e-03 

1.9677 


Example 7.2 Consider the problem (I2.1I) - (I2.3I) in the square domain 17 = (0,1)^. The HWG finite 
element space fc = 1 is employed in the numerical discretization. It has the analytic solution 

^ / -2xy{x - l){y - l)x{x - l){2y “ 1) \ 

\ 2xy{x - l){y - l)y{y - l)(2x - 1) ) 

p = x'^ + y'^ - 

5 


and 
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The right hand side function f in (EU) is computed to match the exact solution. The mesh size is denoted 
by h. 

The numerical results are presented in Tables 7.3-7.4 which confirm the theory developed in previous 
sections. 


Table 7.3. Numerical errors and orders for u of Example 7.2. 


h 

llle.i 

order 

lle.ll 

order 

1/4 

2.8805e-01 


4.2555e-02 


1/8 

1.4913e-01 

0.9498 

1.0184e-02 

2.0630 

1/16 

7.5883e-02 

0.9747 

2.5894e-03 

1.9756 

1/32 

3.8233e-02 

0.9890 

6.5809e-04 

1.9762 

1/64 

1.9173e-02 

0.9957 

1.6598e-04 

1.9872 

1/128 

9.5963e-03 

0.9985 

4.1663e-05 

1.9942 


Table 7.4. Numerical errors and orders for p and A of Example 7.2. 


h 

lk/.ll 

order 

||QbA-A„|| 

order 

1/4 

7.7802e-02 


1.8028e-01 


1/8 

3.7184e-02 

1.0651 

8.2876e-02 

1.1212 

1/16 

1.4725e-02 

1.3364 

3.2275e-02 

1.3605 

1/32 

5.1629e-03 

1.5121 

1.1179e-02 

1.5297 

1/64 

1.6890e-03 

1.6120 

3.5619e-03 

1.6500 

1/128 

5.4636e-04 

1.6282 

1.0720e-03 

1.7324 
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